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ABSTRAGT 

We study a sample of eight massive galaxies that are extreme outliers (3 - 5cr) in the M, - Mbuige 
local scaling relation. Two of these galaxies are confirmed to host extremely large super massive black 
holes (SMBHs), whereas the virial mass estimates for the other six are also consistent with having 
abnormally large SMBHs. From the analysis of their star formation histories and their structural 
properties we find that all these extreme outliers can be considered as relic galaxies from the early 
(z~2) Universe: i.e. they are compact (Re<2 kpc) and have purely old stellar populations (t>10 Gyr). 

In order to explain the nature of such deviations from the local relations, we propose a scenario in 
which the hosts of these uber-massive SMBHs are galaxies that have followed a different evolutionary 
path than the two-phase growth channel assumed for massive galaxies. Once the SMBH and the core 
of the galaxy are formed at 2 , the galaxy skips the second phase, remaining structurally untouched 
and without further mass and size increase. We show that if the outliers had followed the normal 
evolutionary path by growing in size via merger activity, the expected (mild) growth in mass would 
place them closer to the observed local relations. Our results suggest that the SMBH growth epoch 
for the most massive galaxies stopped 10 Gyr ago. 

Subject headings: galaxies: evolution - galaxies: formation - galaxies: stellar content - galaxies: nuclei 
- black hole physics: accretion 


1. INTRODUCTION 


It is now well established that the masses of super 
massive black holes (SMBHs; M, > 10®Mq) strongly 
correlate with some of their h ost galaxies properties, 
such as their luminos i ties ( e.g. iKormendY-^LJiichstoiR 
199,^ iGiiltekin et al.l l200ft iSchulze fc Gebhardtl 120111 


McGonnell fc Ma 20131 iKormendv fc Hd I2013D. their 

masses (e.g Maeorrian et al.l 119981 : iHaring fc Rixl 1^004 

iBeifiori et alT ^12D. or their velocity dispersions (g; 
e.g. ~rrarese fc Merritt! 120001 iGebhardt et al.l 120001 : 


IGiiltekin et al.l I2009I iGraham et al.l 1201111 . This sup¬ 
ports the idea that the host galaxies and their SMBHs 
form and grow in a coordinated way by a com- 


back, mergers or secular evolution (e.g. 

FabianI 119991: 

Honkins et al.1 120061: ISomerville et al.l 12008 

: IMenci et al.1 

20081: l.lahnke & M^ip(jip||20n|l. The correlation between 


SMBH mass and velocity dispersion is extremely tight, 
following a log-normal distribution whose steepness in¬ 
creases as the galaxy sample is enlarged and more accu¬ 
rate measurements of the SMBHs ma sses are used (see 
e.g. th e most recent compila tions of iMcConnell fc Mai 
l 2 ?ril or IKormendv fc Holl2013f l. This correlation implies 
that the most massive galaxies (a '^300-400kms“^) 
should have SMBHs with M, '^IO^^Mq, which in turn 
corresponds to high stellar masses (Mbuige '^ 1 O^^M 0 ) 
and high luminosities (Ly ^ 10 ^^ L 0 ). 

However, in the last few years an increasing num- 
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ber of galaxies have been reported to strongly de¬ 
viate from the above scaling relations. These 
galaxies show larger SMBHs than what it is ex¬ 
pected according to their velocity dispersions or bulge 
masses. This is the case of e.g. the d e eply ana¬ 
lyzed NGG1277 (Ivan den Bosch et al.l [2? )T^ lEmsellepj 


120131 : iFabian et al.l 12013: lYildir im et al J 12015 subm.f l 
NGG4291 and NG G4342 (jBogdan et al.ll2012ll . NGG1332 
(IRusli et al.ll2011h and SPSS ,1151741.75-004217.6 (alias 
bl9. iLasker et all 1201311 . All these extreme outliers lie 
far beyond the intrinsic scatter of the relations, chal¬ 
lenging the supposed co-evolution between SMBHs and 
their host galaxy (>3 (t). The Ho bby-Eberly Telescoye 
Massive G alaxy Survey (hetmgS; Ivan den Bosch etMl 
I 2 OI 2 II 2 OI 5 I) is a new large galaxy compilation aimed at 
searching nearby galaxies in which black hole masses 
can be directly measured. As a by-product, they con¬ 
firmed the existence of two other enormous SMBHs 


in NGG1271 (IWa lsh et a l.l 12015 snbm.f l and MRK1216 
(jYildirim et al.l 12015 subnU i and there ar e few dozens 
more candidates in their sample. Following ILasker et al.l 
(|2013l) notation, these type of extreme SMBHs have been 


named uber-massive SMBHs (UMBHs), although some 
critics have emerged pointing out to an overestimation 
of the real mass. This overestimation could be due to 
the methodology employed (e.g. long-slit spectroscopy 
instead of integral-field t o derive proper dynamical m od- 
els; e.g lEmselleinl 120131 : lYildirim et al.l 12015 subm.!) or 
to variations in the initial mass function (which would 
vary the measured stellar mass of the host galaxy; e.g. 
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lEmsellenil 120131 : iLasker et alJ 120131) . Nonetheless, even 
accounting for these uncertainties, the galaxies hosting 
UMBHs are still extreme outliers of the scaling relations. 
Consequently, it is worth addressing what are the physi¬ 
cal processes that cause these objects to be outliers. One 
possibility is that they are galaxies that did form the 
corresponding stellar mass but then lost i t by e.g. tidal 
stripp ing. However, the recent work from [Bogdan et al.l 
(j2012li . who studied two outliers in the low-mass regime, 
claimed that if that was the case, then the dark mat¬ 
ter halo surrounding the galaxy should also had been 
stripped. Their results, instead, showed that both galax¬ 
ies reside in extended dark matter haloes. 

In this paper we propose and probe an alterna¬ 
tive scenario. We hypothesize that galaxies hosting 
UMBHs form them in an early phase of the galaxy 
evolution (z>2) and once the SMBHs are formed, 
the galaxies remain structurally untouched, which pre¬ 
vents them to reach the present-day scaling rela¬ 
tions, leaving them as outliers. This early forma¬ 
tion is supported by the evolution of th e M, - Mhnip-B 
scaling relation at high redshift (e.e. iWalter et al.l 
2004 l^hnke et ^120091: ICre en e. Peng fc Ludwidl201ClL 
[Petri. Ferrara fc Salvaterra]l2012ll . Therefore, in our sce¬ 
nario, massive galaxies hosting UMBHs are outliers in 
the black hole mass scaling relations because they did 
not follow the expected two-phase growth channel, where 
merger accretion increases the size (and to less extent, 
the mass) of the host gal axies after the central mas¬ 
sive c ore is formed (e.g . iNaab. Johansson fc Ostriked 
120091: lOser et ^ I2010t iHilz. Naab fc Ostriked l2013h . 
Therefore, we would expect the host of UMBHs to 
be galaxies showing the structural properties of the 
massive galaxy population at z^2. In other words, 
they should be the “relics” of the early Universe, 
with old (>10 Gyr) stellar populations and com¬ 
pact (Re <2 kpc) sizes, simil a r to those found at 
high-z (e.g. iTruii l lo et al] l2007t iBuitrago et al.l l2008t 
Carrasco. Conseiicc fc Truiilldl2010l) . The existence of a 
small number of such galaxies without merging accretion 
since z~2 is a natural prediction of the ACDM model (i.e. 
Quilis & Trujillo 2013), but to date only one ca ndidate 
has been robustly confirmed (|Truiillo et al.l[20l^ . 

To probe such prediction, we explore the star formation 
histories (SFHs) for a sample of galaxies with UMBHs. 
If this scenario is realistic, we should be able to pose a 
lower limit for the age of the SMBH formation, which 
corresponds to the age of the host galaxy. The layout 
of this paper is the following. Section 2 describes the 
data used from the different compilations and our spec¬ 
troscopic sample. Section 3 analyses the derived SFHs 
of the sample in order to determine whether a galaxy 
is a relic candidate or not. Section 4 shows the results 
and discusses their implications on the theories of the co¬ 
evolution (or not) of the SMBHs and the host galaxy. A 
final summary is presented in Section 5. Throughout this 
work we adopt a standard cosmological model with the 
following parameters: Hq=70 km s“^ Mpc“^, 17^=0.3 
and Ha=0.7. 

2. THE DATA 

We first collect the largest possible sample of galax¬ 
ies with SMBH mass measurements. There are three 


mai n recently updated com pilations of nearby galax¬ 
ies: iMcConnell fc Mai (I2013L McMlSb iKormendv fc Hoi 
(120131 KH131 and IGraham fc ScottI (1201.4 GSldb They 
are al l based on the initial sample of iGiiltekin et aH 
(l2009fl , but also contain new galaxies and improved mea¬ 
surements on the SMBH masses, galaxy distances and 
bulge mass determinations. All McM13 and GS13 galax¬ 
ies (except for 9 new ones in the latter) are included in 
KH13. In addition, KH13 presents 13 new galaxies. This 
literature compilation corresponds to a total of 95 galax¬ 
ies. 

A number of caveats are worth mentioning. First, 
galaxies found in common among the three compila¬ 
tions (McM13, KH13 and GS13) do not always have the 
same reported values for some of the properties, as dif¬ 
ferent techniques were employed. Second, galaxies are 
morphologically classified following different criteria in 
each compilation. McM13 discriminate between ellipti¬ 
cal galaxies and bulges but do not differentiate between 
classical and pseudobulges, while KH13 make this last 
separation on the bulges, and GS13 only separate by el¬ 
liptical and disk galaxies. In addition, the authors of 
each compilation do not include the same galaxies in 
their fits: some galaxies excluded in a fit may be in¬ 
cluded in another one. Gonsequently, we use here the 
entire compilation to illustrate the correlations (Figure 
1), but we are not deriving new fits for the correla¬ 
tions. In this exercise we use the fits corresponding to 
McM13 for all the galaxies with dynamical mass esti¬ 
mates, with a= 8.32 ±0.05 and /3= 5.64 ±0.32 (for the 
M.-ct; Fig. la) and a= 8.46 ±0.08 and /3= 1.05 ±0.11 
(for M, - Mbuige; Fig. lb). 

Together with the above compilation, we have ex¬ 
panded our sample using the new compilation of galaxies 
from the hetmgs survey. This survey was especially de¬ 
signed to find nearby galaxies that were suitable for black 
hole mass measurements. It is composed of 1022 galax¬ 
ies, hence is an ideal place where to find galaxies with 
very large SMBHs, that will be the basis of our analysis 
for UMBH candidates, as explained in the next section. 

2.1. Defining the UMBH candidates 

To identify UMBH host candidates, we select galax¬ 
ies whose SMBHs are significantly larger than what is 
expected according to the mass of the bulge of its host 
galaxy. In other words, we identify those galaxies which 
strongly deviate from the M, - Mbuige scaling relation. 
We use the following operative definition to select UMBH 
host candidates: a UMBH host candidate is a galaxy 
where the mass of its SMBH is located at a distance of 3cr 
above the M, - Mbuige scaling relation given by McM13 
(see Fig. [T]). With this definition, it is very unlikely that 
the mass of its SMBH can be explained as simply pro¬ 
duced by the uncertainty in the measurement of the mass 
of its SMBH. To have a significant number of UMBH host 
candidates, we use the hetmgs sample. In this paper, 
we are interested in the most massive galaxies, hence we 
will concentrate our efforts on finding UMBH candidates 
in the upper right corner of the M, - Mbuige scaling rela¬ 
tion. 

From the HETMGS sample we first select those galaxies 
that are good candidates to host UMBHs by imposing 
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the following criteria: 

• To have a large dynamical enclosed mass inside the 
central R = 3” aperture of the hetmgs, which 
indicates a very high central mass concentration 
(from either stars or a black hole). We use the 
virial mass estimator Myir = a^RG~^ » 10 ®'^ 
Mq, where G is the gravitational constant and Uc 
the central velocity dispersion. This mass does not 
discriminate between baryonic mass and black hole 
mass. And as such, these virial masses are consid¬ 
ered upper limits for the black hole masses (see 
below). 

• To be nearby enough so that the putative black 
hole can be resolved with high spatial resolution 
facilities; i.e the sphere of influence 9i = G cr^'^ D~^ 
(being D the Hubble flow distance) to b e larger 
then 0.06” 1 see Ivan den Bosch et al.ll2015h . 

With the above criteria, the sample of candidates from 
the HETMGS consists of 173 galaxies. These galaxies 
are, as expected, average massive galaxies with sizes of 
^4kpc. The vast majority of them do not qualify as 
UMBH host candidates (i.e. they are not 3 (t outliers 
of the M,-Mbuige scaling relation). With the 3 (t cri¬ 
teria, our sample is reduced to 30 UMBH candidates. 
Interestingly, this last criteria changes dramatically the 
structural properties of the selected galaxies: they have 
now an average size of ^1.5 kpc. In other words, the 
UMBH host candidates are massive compact galaxies. 
The sizes of the HETMGS galaxies have been retrieved 
from the SDSS database as the half-light effective radii 
in the r-band (re) and then circularized (Re). 

We need to make a last selecting criteria, which is to 
have pre-existing SDS^ spectroscopy with high quality 
(S/N> 20). This is crucial to robustly derive the SFHs of 
the galaxies with the full-spectral-fitting technique em¬ 
ployed in Section 3. This reduces our sample of UMBH 
candidates to 10 galaxies. Two of them (NGC0426 and 
NGC2522) were rejected as their spectra were classified 
by SDSS as Active Galaxy Nucleus (AGN) breadline. 
Having an AGN would not only make the determina¬ 
tion of the stellar population properties of the galaxies 
uncertain but also prevent a proper estimation of the 
structural properties of the galaxies. In fact, the pres¬ 
ence of a bright nucleus could bias the measurement of 
the effective radius of the galaxies towards smaller sizes. 
In order to avoid this potential source of contamination 
from our sample, we have ended with a final sample of 8 
UMBH candidates. 

We use the central enclosed virial mass estimate from 
the HETMGS as an estimation of the M, of the host galax¬ 
ies. This means that both the black hole masses and 
their locus in the SMBH-host galaxy planes should be 
considered as an upper limit. Efforts towards obtain¬ 
ing more accurate estimates using orbit-based dynami¬ 
cal models are underway. For example, the actual black 
hole mass for NGC1277 has been lowered to 12 x 10® 
Mq (jYildirim et al.l 12015 subm.f l , which is a factor of 
1.5 lower than the upper limits. Therefore, if we ap¬ 
ply this factor to the rest of the UMBH candidates we 

^ Sloan Digital Sky Survey 


have a lower, more conservative, estimate of their black 
hole mass. Nonetheless, in the absence of high-resolution 
spectroscopy and orbit-based models for our sample of 
candidates, we will use the upper limit estimates for our 
analysis. However, it is worth noticing that, even consid¬ 
ering the lower estimates, the UMBH candidates are still 
extreme outliers in the local M, - Mbuige scaling relation 
(Fig. lb; ~3-5(T deviations). It is worth noting that, 
when placed in the M, - tr scaling relation, our UMBH 
host candidates are not that extreme (Fig. la; ^ 1 - 3cr 
deviations). The reason for such differences between the 
two planes is further discussed in Sect. 4. The details of 
our complete spectroscopic sample are described in Table 
I. 

2.2. The control sample 

In order to test our scenario, we also explore the ages 
of the stellar populations of galaxies in a control sam¬ 
ple. They were selected from the above (McMI3, KHI3, 
GS13) compilations uniquely by the criteria of having a 
SMBH detection and having high-quality, flux calibrated 
spectra available. No selection based upon morphology, 
size or any other parameter was done. Only 3 galaxies 
out of the 95 have good quality SDSS spectra (NGG3842, 
NGC4261 and NGC4889). High quality archive spec¬ 
troscopic data is also included for a few other control 
galaxi es: NGC3379, NGC4 472. NGC4473 and NGC4697 
(from lYamada et al.l 120061 the reader is referred to this 
publication for a further description of the data). Their 
black hole masses, velocity dispersions and bulge masses 
are those quoted in Table A1 from M cM13. The cir¬ 
culari zed sizes have been retrieved from iKrainovic et al.l 
(|20I3fl f or those galaxies belon ging to the ATLAS3D 
survey (IGannella'ri et all [20Tl]: NGG3379, NGG426I, 
NGC4472, NGG4473 and NGG4697) and from Hyper- 
ledc0 for NGC3842 and NGC4889. This control sample 
is thus composed of large ellipticals and SOs mostly from 
Goma and Virgo clusters. They cover a larger parame¬ 
ter space than our UMBH candidates and are typically 
within the intrinsic scatter of the local relations, showing 
smaller deviations. 

3. ANALYSIS 

As mentioned before, under the proposed scenario 
present-day massive galaxies hosting UMBHs should 
have the characteristics of the typical population of 
galaxies at z^2 after passively evolving since that epoch. 
In other words, they should be massive (M*> 10^^ Mq), 
compact (Re <2 kpc) and with old stellar populations 
(>10 Gyr). Galaxies that fulfill these criteria in the 
local Universe are named relic galaxies and are known 
to b e extreme outliers f rom the local mass-size rela¬ 
tion (jTruiillo et al.l[20?)^ . Therefore, we have quanti¬ 
fied the compactness of our galaxies by measuring how 
much their size deviates from the expected value ac- 
cording to the p resent-day local size-mass relation of 
iShen et ^ (1200311 . This well-known relation was derived 
using a complete sample of 140,000 galaxies from the 
SDS S. The stellar mas s es use d in that paper were taken 
from iKauffmann et al.l (|2003| 1. calculated by multiplying 
the SDSS Petrosian luminosity with the model-derived 

^ http://leda.univ-lyonl.fr 
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TABLE 1 

The spectroscopic sample 


Galaxy 

Spectra 


Mbulge( + i —) 

vel. disp. 

Re 

deviation (+, —) 

deviation (+, —) 



(10® Mq) 

(IQii Mq) 

(kms“^) 

(kpc) 

^bulee 

vel. disp. 

UMBH host candidates from HETMGS 

NGC1270 

SDSS 

<12 

3.3 (2.7,1.3) 

353±4 

1.92±0.03 

3.2a (-,1.6) 

1.0(7 (-,0.7) 

NGC1271 

SDSS 

<8 

1.5 (0.8,0.6) 

293±6 

1.31±0.02 

4.0a (-,2.0) 

2.0(7 (-,1.9) 

NGC1277 

SDSS 

<17 

2.1 (1.9,0.7) 

385±4 

1.64±0.03 

4.2a (-,1.4) 

0.8(7 (-,0.7) 

NGC1281 

SDSS 

<5 

1.2 (1.0,0.4) 

258±4 

1.46±0.03 

3.4o- (-,0.6) 

2.0(7 (-,0.9) 

NGC2767 

SDSS 

<5 

1.4 (1.1,0.6) 

239±3 

1.95±0.03 

3.2o- (-,1.7) 

2.5(7 Li-3) 

PGC012557 

SDSS 

<7 

0.6 (0.6,0.3) 

283±4 

0.71±0.02 

4.7a (-,1.0) 

1.8(7 (-,0.8) 

PGC012562 

SDSS 

<6 

0.4 (0.3,0.2) 

259±4 

0.70±0.02 

5.0(7 (-,1.2) 

2.2(7 (-,0.7) 

PGC032873 

SDSS 

<14 

2.3 (1.2,0.9) 

356±3 

1.76±0.03 

3.7a (-,1.2) 

1.0(7 (-,0.8) 

Control sample from McM13 

NGC3379 

Y06 

0.4 (0.1,0.1) 

0.7 (0.5,0.3) 

206± 10 

3.3±0.4 

1.0(7 (1.1,1.1) 

0.6(7 (0.6,0.8) 

NGC3842 

SDSS 

9.7 (3.0,2.5) 

16.0 (14.0,7.0) 

270± 14 

14.0±1.3 

0.8(7 (1.1,1.2) 

2.4(7 (0.7,0.9) 

NGC4261 

SDSS 

0.5 (0.1,0.1) 

8.3 (7.0,4.0) 

315±15 

7.4±0.7 

2.1(7 (1.1,1.1) 

1.9(7 (0.5,0.8) 

NGC4472 

Y06 

2.5 (0.6,0.1) 

9.0 (7.0,4.0) 

300± 15 

8.2±0.4 

0.2(7 (0.8,1.1) 

0.2(7 (0.6,0.7) 

NGC4473 

Y06 

0.1 (0.0,0.0) 

1.6 (2.0,0.6) 

190±9 

4.3±0.3 

2.1(7 (0.9,1.6) 

0.6(7 (0.7,1.0) 

NGC4697 

Y06 

0.2 (0.0,0.0) 

1.3 (0.7,0.6) 

177±8 

6.4±0.4 

0.8(7 (1.3,1.5) 

0.7(7 (0.8,1.5) 

NGC4889 

SDSS 

21.0 (16.0,15.5) 

18.0 (14,8) 

347± 17 

12.6±1.4 

1.6(7 (1.6,2.5) 

1.7(7 (1.0,2.6) 


Column 1 - Our spectroscopic sample from the compilation of hetmgs (top) and McM13 (bottom). Column 2 - origin of the spectra: 
SDSS or Y06 (from Yamada et al. 2006). Columns 3 and 4 — the masses of the SMBHs and the masses of the bulge of the galaxy from 
their dynamical estimates (from HETMGS and McM13). Column 5 - Velocity dispersions (derived with STARLIGHT and from McM13). 
Column 6 — Circularized effective radii (from the SDSS r-band, ATLAS3D and Hyperleda). Columns 7 and 8 - Deviations from the local 
SMBH mass scaling relations. Note that for the UMBH host candidates the black hole masses are an upper limit, thus marked with an <. 



(7 (km/s) Mbuige/M© 

Fig. 1.— Black hole mass scaling relations with velocity dispersion (left) and as function of bulge mass [right). Filled black diamonds 
correspond to the 8 UMBH host candidates, while circles correspond to the compilation of published galaxies with SMBH detections (grey), 
with our 7 control galaxies emphasized in black large circles. The black solid line is the fit corresponding to McM13 and the dashed lines 
limit the Ic, 3f7 and 5 (t deviations from it. Arrows mark the expected growth for individual massive galaxies since Z'^2 if they had followed 
the merging phase of assembly for massive galaxies, making them more normal in terms of SMBH mass (see section 4). The lower error 
bars show the ~1.5 factor that the black hole mass is estimated to possibly vary using more accurate techniques. 
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Fig. 2. — Stellar mass-size distribution of the SDSS galaxies. 
Overplotted in green, the Shen et al. (2003) mass-size distribution 
for SDSS n >2.5 galaxies is shown. The position of our UMBH 
host candidates is shown with black diamonds, while the control 
sample is shown in blue filled dots. It is clear that all our UMBH 
host candidates are extreme outliers of the local mass-size relation 
and thus good candidates to be massive relics of the early Universe. 

mass-to-light ratio. The radii considered were the Sersic 
half-light radius i? 5 o,s in the z-band. For galaxies with 
n > 2.5, the local mass-size relation is: 

R{kpc) = h{M/MQY (1) 

where R is the mean radius, M is the stellar mass, and 
a and b are fitting parameters given by a least-squares fit 
to the data, a= 0.56 and 5=2.88x 10“®. This relation 
is shown as the solid green line in Figure 2. This fig¬ 
ure indicates the location of our spectroscopic sample in 
the mass-size relation plane. It is clear that our UMBH 
candidates largely deviate from the local relation. 

We consider here that a galaxy is compact if the ra¬ 
tio Re/R.e.shen 03 ^0.33 (see Table 2). This corresponds 
to ^1/3 of the normal size that massive galaxies show 
at z = 0 at a given stellar mass. The location in the 
stellar mass-size plane of our UMBH host candidates is 
equivalent to the locus occupied by the massive compact 
galaxies at high-z. 

To determine the stellar populations of both the out¬ 
liers and the control sample we derive the SFHs of 
our galaxies applying a full-spectral-fitting approach. 
This technique has several advantages compared to the 
classical line strength approach. One of them is that 
the results are not limited to luminosity-weighted es¬ 
timates but also give the output in terms of mass, 
allowing to recover the true fossil imprints of the 
epoch when the stars wer e created. The code em¬ 
ployed here is STARLIGHT (iCid Fernan des et al.l 120051 : 

I Cid Fernandes fc Gonzalez Delgadoll2010ll . which creates 
a combination of single stellar populations (SSPs) model 
predictions that best resemble the observed galaxy spec¬ 
trum while minimizing the The reader is referred 
to the source code papers for a deeper description of 
the methodology, or to the Appendix for a brief sum¬ 
mary of this technique. Th e SSP models used ar e 
MIUSCAT (|Vazdekis et al .11^0121 : iRicciardelli et all 1201211 . 
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which cover a wide range of both ages (0.1-17.8 Gyr) 
and metallicities (-1.71 to 0.22dex). For this exer¬ 
cise in particular, 46 ages and 4 different metallicites 
were considered, making a total of 184 SSP templates. 
These SSP models also allow for variations on the 
IMF slope and shape. The IMF h as been recently 
reported to be non-universal (e.g. iCappellari et al.l 
20121: iConrov fc van Dokkumll20i2l: iSpiniello et al.ll2012l : 


Ferreras et al.ll2013t iLa Barbera et al.ll2013D and to have 
a stro ng impact in the derived stellar populatio n prop¬ 
erties J^T^Mateu^^MdekIF^^pJaTos3[2Qi3). How¬ 
ever, iMartfn-Navarro et al.n~oT^ have shown recently 
that for the compact massive relic galaxy NGC1277, 
changing the IMF does no modify the derived age prop¬ 
erties. This indicates that it is still unclear which is the 
parameter that really drives these IMF variations (e.g. 
velocity dispersion, metallicity, a-enhancement). There¬ 
fore, here we will focus on the standard assumption of a 
Kroupa Universal IMF. 

Figure 3 illustrates the derived SFH from the spectral 
fitting for one of the galaxies (NGC1271). The derived 
SFHs for the rest of our spectroscopic sample can be 
found in the Appendix. The right panels in Figure 3 
show the derived SFHs as the percentage of mass cre¬ 
ated at a given age, considering both the contribution of 
each individual metallicity from the models and of the 
total metallicity. It is worth emphasizing here an im¬ 
portant caveat o n the form the derived SFH s are shown. 
As mentioned bv iGid Fernandes et al.l (|2nn5h . the reader 
should not take the form of each single burst at face 
value. This is partly because the SSP models do not 
cover uniformly the entire age parameter space. For ex¬ 
ample, there is always a gap at 13 Gyr because there is 
no SSP model corresponding to that age in the MIUSCAT 
models. Therefore, the derived SFHs should be better 
interpreted as clearly differentiated epochs of formation, 
i.e. young, intermediate and old. These considerations 
are fundamental for the purpose of our exercise, as one 
of our criteria to consider a galaxy as a relic of the z~2 
Universe is that all their stellar populations must be old 
(>10Gyr). When looking at the age distributions, it 
can be seen that sometimes there is a small contribution 
at 8-9Gyr, which can be considered as the tail of the 
distribution that peaks at ~ 10 Gyr. For this reason, in 
order to be more quantitative in our statement whether 
a galaxy has or not entirely old stellar populations, we 
measure the fraction of mass H which is within the inter¬ 
val 8<t<10 Gyr. We thus consider a galaxy compatible 
with having all its stellar mass formed at z> 2 if there 
are no stellar populations younger than 8 Gyr and H is 
below 10% of the total mass. 

We summarize in Table 2 the quality parameters from 
the fit and whether the outliers satisfy our criteria to 
be considered relic galaxies or not. We find that 7 out 
of the 8 objects undoubtedly qualify as relic galaxies, 
while only one is dubious. As these compact relic galax¬ 
ies are also outliers in the present-day mass-size relation 
(Figure 2), finding these UMBH host candidates could 
constitute a new way of detecting these rather elusive 
relic ga laxies in the nearby U niverse (e.g. iTruiillo et al.l 
20091: PTavlor et’all l2010l: iFerre-Mateu et al.l ^121: 


Truiillo. Garrasco fc Ferre-MateuI 120121 iPoggianti et all 
20131 1. Regarding our control sample, none of the galax- 
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ies can be considered as massive relic of the high-z Uni¬ 
verse because of their large sizes. However, a couple of 
them are as old as our relic galaxies, indicating that at 
least the stellar populations of their most central parts 
were formed at very early epochs too. We would like to 
emphasize that our criteria based on the fl value, in order 
to be fully consistent with the statement that the galaxy 
is a relic, should be shown to be true along the entire 
structure of the galaxy. As we are employing SDSS data, 
the coverage of the fiber is ~0.3-0.8Re for the compact 
galaxies and ^0.01-0.03Re for the control sample. Con¬ 
sequently, the SDSS data does not allow us to perform 
such radial study. However, thi s issue has been asse ssed 
for the compact candidates in iTruiillo et al.l (I2014H for 
NGC1277 and in the companion paper of Ferre-Mateu 
(2015b, in prep.) for PGC032873. These works both 
use high-quality spectroscopy from the William Herschel 
Telescope to the derive SFHs and stellar populations 
properties out to several effective radii, showing no gradi¬ 
ents and supporting our hypothesis of this type of galax¬ 
ies being relics. 


4. DISCUSSION 

The tight correlations that the mass of the SMBHs and 
their host galaxy properties follow have been considered 
as a proof of the co-evolution between galactic bulges and 
the SMBHs they host. According to this paradigm, they 
should be coupled by a common physical mechanism, e.g. 
feedback, mergers, secular evolution. However, a small 
sample of galaxies hosting UMBHs now challenge this as¬ 
sumed universal co-evolution as they clearly fail to follow 
such relations. In this paper we have proposed that the 
nature of these outliers hosting UMBHs is connected to 
the uncommon evolutionary path followed by these mas¬ 
sive galaxies (without growing in mass or size). It is gen¬ 
erally assumed from galaxy formation models that mas¬ 
sive galaxies form in a two-phase mechanism, where the 
central massive part of the galaxy is created in a fast and 
violent event at high-z, and then the galaxy grows by pos¬ 
terior merger activity fe.g. iNa ab. Johansson fc OstrikeJ 
120091 lOser et al.l[2010HHilz. Naab fc Ostrikeil2013ll . But 
because mergers are stochastic events, it is expected that 
some galaxies remain untouched over cosmic time, re- 
maining compact and entirely old (see iQuilis fc Truiill^ 
I2013f) . Therefore, if the UMBH-outliers are those relics, 
as our results suggest, then the nature of the deviations 
could be explained by the lack of merging activity and its 
consequent lack of galaxy growth. Under our hypothesis, 
these galaxies should occupy in the SMBH mass scaling 
relations the position of the population of galaxies at 
z^2. Gonsequently, present day most massive galaxies 
should start at those locations ^10 Gyr ago and their 
growth in mass and size should be able to move them 
to the current scaling relations. If this is true, the ex¬ 
pected mass and a increase of the massive galaxies since 
z^2 should be enough to locate our outliers in agreement 
with the local relations. 

We explore this possibility on what follows. The ar¬ 
rows in Figure 1 show the estimated growth in veloc¬ 
ity dispersion and stellar mass since z^ 2 for individ¬ 
ual massive galaxies. From semi-analytical models, if 
an individual galaxy grows ^7 times in size , it can in- 
crease its mass by almost a factor of 5 (e.g. lOser et ^ 


[Mnl: lTruiiHo7 Ferreras fc de La Rosal[20Tlll . This could 
account for the missing stellar mass, placing the UMBH- 
outliers closer to the central distribution in the M,- 
Mbuige relation. As commented above, our sample is not 
that extreme when placed in the M, - a relation, hence no 
much variation in the velocity dispersion should be ex- 


dictions fe.g. 

Genarro & Truiilloll2009HOser et al.ll2012t 

lOosfi & Habd 

20131 IWellons et al.l 1201511. where the ve- 


locity dispersion varies very slightly since z^ 2. Depend¬ 
ing on the dominant merger channel, individual galaxies 
can increase very mildly thei r velocity dispersion by a fac¬ 
tor of ^1.1 (|Hilz et al.ll2012l : Tapia et al. 2015 subm.). It 
is worth noting that in our scenario, the arrows move hor¬ 
izontally because an increase in the SMBH is expected 
to be almost negligible during the second phase. It is 
typically assumed that the amount of merge rs that oc¬ 
cur s i nce z^ 2 have a mass ratio of 1:3 to 1:5 (|Oser et al.l 
120101 [ 201 ^ . Here we assume that the massive galaxies, 
once formed at z^2, have a much larger SMBH than 
what they should have from their stellar masses. Then, 
they merge with smaller satellite s, whose mass-size e volu- 
tion is milder over cosmic time (|Truiillo et al.l[2?)?)3 l and 
therefore are expected to be located already in, or very 
close to, the local relation having smaller black holes. 
This type of encounters would produce a negligible in¬ 
crease in the mass of the SMBH, as se en in cosmologi- 
cal simulations of S MBH growth (e.g. lYoo et al]l2007t 
iWellons et al.l[20T^ iKnlier et al.l 1201 51 1. Therefore, it is 
a valid assumption to account for evolution only in the 
Mbuige and a directions for the relic galaxies and neglect 
the one in the SMBH mass direction. As a final comment, 
considering the lower black hole estimates would further 
support our scenario: the UMBH candidates would still 
be extreme outliers in the M,-Mbuige relation but they 
would be located even closer to the local scaling relations 
after accounting for the missing mass. 

In a nutshell, our scenario is schematically depicted in 
Figured. The left panel shows the normal evolution that 
a massive galaxy follows over cosmic time to become a 
large ell iptical galaxy to day, with the two-phase of for¬ 
mation (|Oser et al.ll201Clf) . The right panel illustrates the 
expected evolutionary track for the relic/outlier galaxies. 
We remind the reader that this is a cartoon to illustrate 
the proposed scenario, where a number of assumptions 
have been here simplified, as many questions in the (co- 
) evolution of SMBHs and their host galaxies remain un¬ 
solved. For example, the exact position in the early Uni¬ 
verse phase. It is not yet clear if the SMBH and the 
host started forming at the same time, but several ob¬ 
servational and theoretical works suggest that the SMBH 
started forming first, as inferred from the mild deviations 


Jahnke et al.l 120091 IMerloni et al.l 120101 

Heines et al. 

2011 

; IPetri. Ferrara & Salvaterral 120121: 

Khandai et al. 


2 OI 2 I I. Our scenario also seems to support this case, as 


by the end of the z> 2 phase the SMBH should be almost 
fully in place, being larger than expected from the galaxy 
stellar mass. Although we are not in the position to con¬ 
firm this from our derived SFHs, we can pose a lower 
limit for the SMBHs growth at ^ 10 Gyr, when the star 
formation activity in the host galaxy is halted. At this 
point, one would expect the massive compact galaxy to 
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Fig. 3. — Illustration of the full-spectral-fitting approach used in this work for one of our UMBH host candidates. The whole sample can 
be found in the Appendix. The left panel shows the galaxy NGC1271 SDSS spectrum (black), the fitting from STARLIGHT (green), and the 
residuals from the fitting and the masked spectral regions to avoid possible emission line contributions (grey). The right panel shows the 
mass-weighted derived SFHs considering both the contribution from each individual metallicity and from the total metallicity. 


TABLE 2 

Fitting and Massive galaxy relic criteria 


Galaxy 

[x1 

adev (%) 

R'e/Re.shenOS 

%M(t < 8Gyr) 

n{%) 

Relic? 

UMBH host candidates 

NGC1270 

1.6 

1.7 

0.24 

0 

7 

/ 

NGC1271 

1.7 

1.6 

0.25 

0 

10 

/ 

NGC1277 

1.9 

1.5 

0.26 

0 

0 

/ 

NGC1281 

1.2 

1.6 

0.33 

0 

0 

/ 

NGC2767 

0.7 

1.8 

0.37 

0 

0 

/ 

PGC012557 

1.3 

1.7 

0.22 

0 

0 

/ 

PGC012562 

1.3 

1.3 

0.23 

0 

14 


PGC032873 

1.0 

1.0 

0.23 

0 

0 

/ 

Control 

NGC3379 

0.1 

1.9 

0.92 

5 

4 

7 

NGC3842 

0.8 

2.0 

0.73 

1 

14 

X 

NGC4261 

0.7 

1.8 

0.55 

0 

5 

X 

NGC4472 

1.5 

1.7 

0.55 

7 

5 

X 

NGC4473 

1.1 

1.3 

0.83 

1 

14 

X 

NGC4697 

1.2 

1.4 

1.26 

21 

7 

X 

NGC4889 

0.8 

1.8 

0.61 

0 

10 

X 


The table summarizes the different criteria by which a galaxy can be considered a relic candidate. Column 2 and 3 are measurements of 
the quality of the fit from STARLIGHT, (see Appendix). Col umn 4 measures th e compactness of a galaxy, by showing how much the galaxy 
size deviates from the present-day mass-size relation of IShen et al.l II2003I1 for n > 2.5 galaxies. Column 5 and Column 6 quote the 
fraction of stellar populations with t < 8Cyr, and the fraction of mass with 8<t<10 Cyr, f2. Column 7 states the verdict on whether the 

galaxy is considered a relic candidate or not. 


move into the second phase of formation, balancing out 
the deviations with its SMBH, which roughly evolves, by 
growing largely in size and mildly in mass through merg¬ 
ers. This will eventually place the massive galaxy in the 
local M, - Mbuige relation. 

We are aware of the limited sample employed here, 
but so is the published sample of galaxies with de¬ 
tected SMBHs and accurate measurements or the sam- 
ple of relic galaxy candidates. In a recent effort, 
ISaulder. van don Bosch fc Mieskel ()2015D have revisited 
the number of relic candidates i n the SDSS, enlarging 
the one from iTruiillo et all (|2009D to 76 candidates. The 
next step is to find out how many of these galaxies host 
a SMBH and if so, if they are UMBHs. This will be 
achieved with the upcoming era of the largest telescopes, 
when SMBH detections and measurements will be im¬ 
proved. On the meantime, surveys such as the het- 
MGS are ideal compilations for new SMBH detections. 
A deeper understanding on the formation of these ex¬ 
treme outliers (both in the upper and the low-mass end) 


is crucial to test the current theories of galaxy formation 
and evolution and cosmological models. 

5. SUMMARY 

In this work we have studied a sample of 8 compact 
massive galaxies selected from the hetmgs that are can¬ 
didates for hosting UMBHs and that are therefore con¬ 
sidered extreme outliers from the correlation that the 
SMBHs and the mass of their host galaxies follow. We 
propose that they deviate because the host galaxies did 
not structurally evolve in the way expected for massive 
galaxies due to their relic nature. Our main results are 
here briefly summarized: 

• We find that 7 out of the 8 UMBH host candidates 
fulfill the criteria to be considered a relic galaxy 
(i.e. Re<2 kpc and t>10 Gyr), while one of them 
is a dubious case. Therefore, selecting galaxies that 
are extreme outliers in the SMBH mass scaling re¬ 
lations could represent a new way to find the elusive 































A. Ferre-Mateu et al. 


I09 Core of Phase II: nassive galaxy M 

(SMBH mass) massive galaxy grows via mergers (SMBH mass) 


/ 


O •• 


Today 


z~2 




W ^ . 




l5il 


Early Universe 


Massive galaxy evolutionary track 


mass) 


Core of 

massive galaxy Mo Phase II: 

is formed no merger adWty and 
stnctaral evolution 


/ 

> 

•s.* C 

/ 



z~2 anjl^^oday 


log 

-nCalaxy 


Early Universe 


ReUc galaxy evolutionary track 


log 

’ (Galaxy 
mass) 


Fig. 4. — Schematic view of the scenario proposed in this work: galaxies hosting UMBHs depart from the mass scaling relationship due 
to their unusual growth channel. The figure shows the assumed evolutionary track over cosmic time in the SMBH mass- host galaxy mass 
plane for massive galaxies following the two-phase growth channel (left) and for the relic galaxies (right)^ who skip the second phase of 
such channel of, formation, Tremaining unaltered over cosmic time. ^tt i t ,i mv /ttvtt 

relic galaxies in tne nearoy Universe. • We can pose a lower limit to the age on the SMBH 

growth at ^ 10 Gyr, based on the derived SFHs. 


• When we plot the M.-cr and M,-Mbuige relations, 
the loci of the galaxies hosting UMBHs should rep¬ 
resent the position of the massive galaxy popula¬ 
tion at z^ 2 . As cosmic time progresses, galaxies 
following normal growth paths should evolve to¬ 
wards galaxies with larger host masses in those 
planes. 


This work was supported by the Japan Society for 
the Promotion of Science (JSPS) Grant-in-Aid for Sci¬ 
entific Research (KAKENHI) Number 23224005 and by 
the Spanish Ministerio de Economia y Gompetitividad 
(MINEGO; grants AYA2009-11137, AYA2011-25527 and 
AYA2013-48226-C3-1-P). 


APPENDIX 

STARLIGHT AND THE STAR FORMATION HISTORIES 

We briefly summarize here how the full-spectral-fi tting code STARLIGHT works, but the reader is referred to the 
sourc e code papers for a more detailed explanation (|Cid Fernandes et al.l I2005L iGid Fernandes fc Gonzalez Delgadol 
I2010H . The code models the extinction as due to foreground dust, and different reddening-laws can be selected to 
correct from Galactic extinction. Then, it finds the fraction Xj that a given jth SSP contributes to the total flux of the 
galaxy (normalized to a certain wavelength Aq), creating a synthetic spectrum. In other words, it creates a combination 
of SSPs that best resemble the observed galaxy spectrum and that minimize the The first two parameters presented 
in Table 2 are indicative of the quality of the fit. [x^] represents the total x^ divided by the number of A’s used in the 
fit. The second parameter, adev, is a proxy for the mean deviation over all fitted pixels. Good values for a fit are those 
with adev below 2-3%. A visual inspection of the fits and their residuals confirms the quality of the results obtained 
(see Figures Al and A2). 

From each Xj we can derive its contribution in mass, 7 ^, directly from the mass-to-light ratio of each individual SSP 
(from the models), to obtain the desired mass-weighted estimates such as the star formation episodes and mean stellar 
populations parameters, which are derived as: 


(G)m = tj ; for the age (Al) 

{Z.) M = E jjZj ; for the metallicity (A 2 ) 

1=1 


We can therefore reconstruct the SFH of the galaxy by summing up all the fractions of stellar mass created at a given 
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epoch, as shown in Figures A1 and A2. They present the whole spectroscopic sample data and their derived SFHs. 
Like in Figure 3, the left panels show the galaxy spectrum with fit from STARLIGHT and the right side ones, the derived 
SFHs in terms of stellar mass. 
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Fig. 5.— Star formation histories derived for the sample of UMBH host candidates. Each row corresponds to a galaxy. Left panel 
shows the galaxy spectrum (black), the fitting (green), and the residuals from the fitting and the masked spectral regions to avoid possible 
emission line contributions (grey). The right panels show the mass-weighted SFHs considering both the contribution from each individual 
metallicity and from the total metallicity. 
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Fig. 6. 
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— As in Figure Al but now for the control sample. 














































